Coexistence of Sharp Quasiparticle Dispersions and Disorder Features in Graphite 
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Angle resolved photo-electron spectroscopy (ARPES) on azimuthally disordered graphite demon- 
strates that sharp quasiparticle dispersions along the radial direction can coexist with a complete 
lack of dispersion along the azimuthal direction. This paradoxical coexistence can be explained in 
terms of van Hove singularities in the angular density of states. In addition, non-dispersive fea- 
tures at the energies of band maxima and saddle points are observed and possible explanations are 
discussed. This work opens a new possibility of studying the electronic structure of novel layered 
materials using ARPES even when large single crystals are difficult to obtain. 
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The ability to sharply resolve crystal momentum val- 
ues of single particle excitations has made ARPES a very 
powerful tool in addressing the electronic structure of 
solid, as has been successfully demonstrated on single 
crystalline samples over the past decades jl| . Due to the 
translational symmetry along the surface of a single crys- 
tal, the crystal momentum parallel to the surface (k|| ) is 
conserved during the photoemission process, allowing a 
complete momentum space map of the initial state. This 
holds despite the short photoelectron lifetime 0, 0, H 
which can severely broaden the resolution of the momen- 
tum perpendicular to the surface (k z ). Indeed, even in 
the limit of an extreme k z broadening that results in no 
resolution of k z , strong ARPES dispersions are expected 
as a function of k| , since the one dimensional density of 
states (ID-DOS) D Z (E) oc dk z /dE obtained by integrat- 
ing over k z is dominated by contributions from van Hove 
singularities in high symmetry planes [4| . 

On the contrary, for those systems characterized by 
orientationally disordered domains, i.e. polycrystalline 
materials, the translational symmetry is preserved only 
within each domain. As a consequence, the dispersion 
measured by ARPES is the average dispersion over dif- 
ferent domains, or equivalently azimuthal angle (f>, which 
in general leads to no dispersion. However, extending 
the ID-DOS D Z (E) scenario for k z discussed above fur- 
ther to the plane, there is an interesting possibility that 
a layered polycrystalline sample, with a strong azimuthal 
disorder, can nevertheless give distinct dispersions in the 
radial direction. This would happen if the average dis- 
persion is dominated by those along the high symmetry 
directions due to van Hove singularities in the angular 
density of states D^(E) oc d0/dE. To date, this possibil- 
ity has never been demonstrated experimentally and pho- 
toemission studies on disordered samples have focused on 
angle-integrated features without any momentum infor- 
mation. 

In this paper, we report a high resolution ARPES 
study on the electronic structure of azimuthally disor- 
dered graphite. For the hrst time, we report clear evi- 
dence that sharp quasiparticle dispersions, in agreement 
with band structure calculation along the high symme- 



try directions, can coexist with a circular Fermi energy 
intensity map, a definitive signature of azimuthal disor- 
der [5|. In addition we report non-dispersive features 
at the energies of band maxima or saddle points, which 
are attributed to the loss of momentum information by 
indirect photoemission process or elastic scattering. A 
practical implication of this study is that more ARPES 
opportunities can be made available for layered materi- 
als even when high quality single crystals of large size are 
difficult to obtain. 

ARPES data were collected at beam line 10.0.1 of the 
Advanced Light Source (ALS) at the Lawrence Berkeley 
National Laboratory, using an SES-R4000 analyzer. The 
wide angular mode with acceptance angle of 30° and an- 
gular resolution of 0.9° was utilized for most scans, while 
high resolution angular mode with acceptance angle of 
14° and angular resolution of 0.1° was utilized for one 
scan. The total instrumental energy resolution was 15 
meV at 25 eV photon energy and 25 meV for other pho- 
ton energies used (40, 55, 60 eV). The sample used was a 
grade ZYA highly oriented pyrolytic graphite (HOPG), 
obtained commercially from Structure Probe Inc. The 
sample was cleaved in situ in an ultra high vacuum bet- 
ter than l.OxlO -10 Torr and measured at temperature 
50 K. 

Figure 1(a) shows an ARPES intensity map mea- 
sured at the Fermi energy taken at 40 eV photon en- 
ergy. Throughout this paper, we use a color scale such 
that black represents high intensity in the raw data and 
blue represents low intensity. According to band struc- 
ture calculation, a constant k z cross section of graphite 
Fermi surface can be a small hole pocket, a small elec- 
tron pocket, or a point located at the six corners of the 
hexagonal Brillouin zone (dashed line in figure 1(a)), de- 
pending on the value of k z 0, 0j- Experimentally, the 
predicted small electron or hole pockets have yet to be 
resolved, and measurements on single crystalline samples 
have shown only small dots of high intensity at these cor- 
ners |l| , schematically drawn as shaded circles in Figure 
1(a). For the graphite sample under study, the Fermi 
energy intensity map, symmetrized by three fold rota- 
tions to fill the entire Brillouin zone, shows a perfectly 
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FIG. 1: (a) Fermi energy intensity map. The hexagonal Brillouin zone (dashed lines) and Fermi surface (shaded circles) expected 
for single crystalline graphite are drawn schematically, (b) Intensity map versus binding energy and in-plane momentum along 
the solid line in (a) taken at 60 eV photon energy. Arrow marks the Fermi energy crossing point (kf). The inset shows EDC 
at ki? taken at 25 eV photon energy in the high angular resolution (0.1°) mode, (c) Second derivative of raw data in (b) with 
respect to energy. LDA dispersions along both F-K-M' direction (solid lines) and F-M-r' direction (dashed lines) are plotted 
for comparison. The Brillouin zones are labeled on top of this panel for the two high symmetry directions. 



circular pattern |18j . in contrast to what is expected for 
single crystalline graphite. This is attributed to the an- 
gular spread of the dots to a circle due to the azimuthal 
disorder of the sample 

Figure 1(b) shows an ARPES intensity map as a func- 
tion of binding energy and in-plane momentum k|| , cor- 
responding to the momentum cut shown as a solid line in 
Figure 1(a). Despite the strong azimuthal disorder giv- 
ing a circular Fermi energy intensity map in Figure 1(a), 
we observe, surprisingly, very clear dispersions over the 
entire energy range. Furthermore, at the Fermi energy 
crossing point kp, a sharp coherent quasi-particle peak 
is observed. This is shown in the inset, where an en- 
ergy distribution curve (EDC), energy cut at a constant 
momentum, is plotted. Here the half width of the EDC 
peak is 20 meV (50 meV FWHM due to the asymmetry 
of the line shape), defining the sharpest peak observed in 
graphite so far Q. 

In Figure 1(c) we report the second derivative of the 
raw data of Figure 1(b) with respect to energy. The sec- 
ond derivative method has been used in the literature to 
enhance the direct view of the ARPES dispersion. Lo- 
cal density approximation (LDA) band dispersions along 
two high symmetry directions T-K-M (solid lines) and 
F-M-F (dashed lines) are plotted in the same figure for 
a direct comparison. Despite a polycrystal-like sample 
implied by the Fermi energy intensity map, an excellent 
agreement is observed between the experiment and the 
theory. We can identify the dispersions between 4 eV and 
23 eV as originating from the sp 2 orbitals with strong 
intra-layer a bonding (black lines), and the dispersions 
between Fermi energy and 11 eV as originating from the 
p z orbitals with weaker 7r bonding (white lines). We note 
that the calculated dispersions were stretched by 20% in 
energy throughout this paper, as suggested in the litera- 



ture 0, |U 0| . The stretching of the LDA band disper- 
sions is attributed to missing self-energy corrections in 
LDA, since ab initio quasiparticle calculations based on 
the GW method show that for graphite the quasiparticle 
band dispersion near the Fermi level is 15% larger 9]. 

The direct comparison between Figure 1(a) and Figure 
l(b,c) shows an apparent paradox in our data, namely, 
the coexistence of azimuthal disorder feature (Figure 
1(a)) with single crystalline features (Figure l(b,c)). This 
can be readily understood if we consider an angular av- 
erage of the calculated dispersions. Such an angular av- 
erage would be necessary if the sample consisted of many 
small single crystallites with strong azimuthal disorder. 

Figures 2(a-c) show ARPES cuts for three azimuthal 
angles, 0=0°, 10°, 20°. A direct comparison between 
panels a, b and c shows no appreciable angular de- 
pendence of the dispersions, establishing that the az- 
imuthally invariant electronic structure of Figure 1(a) 
at the Fermi energy extends to the entire band width. 
Thus, these data strongly support the azimuthal disor- 
der model described in the previous paragraph, and in- 
dicate that the ARPES data measured are actually a 
ID-DOS D^(E) along the azimuthal direction <fi, in anal- 
ogy with the well-known ID-DOS D Z (E) along the k z di- 
rection Q . As in the latter case, then, one would ex- 
pect that van Hove singularities arising from states along 
the high symmetry directions to contribute dominantly, 
and this gives an explanation why the measured disper- 
sions accurately reflect the dispersions along the two high 
symmetry directions. The ID-like van Hove singularities 
arise from states along high symmetry lines because these 
states have zero group velocity along the arc with con- 
stant k|| magnitude. 

In Figure 2(d, e) we show LDA calculations support- 
ing this reasoning. In panel d, we show the dispersion 



3 




FIG. 2: Dispersions for azimuthal angles (f>=Q° (panel a); 10° (panel b) and 20° (panel c). The <f> angle is defined in the inset of 
panel a. LDA band dispersions along T-K-M' (solid lines) and r-M-F' (dotted lines) are plotted for comparison, (d) Calculated 
dispersions for single crystalline graphite along an arc (shown in the inset) with radius equal to FK distance, (e) Calculated 
density of states D^(E) for single crystalline graphite by integrating over an arc from A to B (see inset of panel d). Singular 
peaks in D^(E) occur at energies corresponding to band energy extrema, some of which are shown in panel(c,d) as shaded 
circles and open circles for A (FK direction) and B (rM direction) respectively. 



for single crystalline graphite along an arc from A to 
B with radius equal to TK distance. Within this arc, 
the TK direction corresponds to point A and the TM di- 
rection to point B respectively. As expected, extrema 
in the band dispersion occur at the two high symme- 
try directions, A (shaded circles) and B (open circles). 
The calculated ID-DOS D^(E) over this arc, and thus 
over the entire azimuthal angle range by symmetry, is 
shown in panel e. One can see diverging ID van Hove 
singularities as sharp peaks occurring at energies where 
bands cross points A and B, which completely dominate 
over other contributions. This nicely explains why well- 
defined sharp peaks with large dispersions can be ob- 
served in this azimuthally disordered sample, despite the 
fact that the observed data come from averaging over all 
azimuthal directions. 

The data presented so far can be summarized as show- 
ing well-defined dispersions along the radial direction 
with a complete lack of dispersion along the azimuthal 
direction. Therefore, our data suggest that the graphite 
sample under study consists of finite size single crystalline 
grains much smaller than the analysis area (« 100 /im) 
with a complete azimuthal disorder. However, each grain 
is large enough to allow for highly dispersive quasiparti- 
cles to exist. In addition, we have measured the disper- 
sion perpendicular to the surface, k z dispersion, using 
photon energy range from 34 to 155 eV at beam line 
12.0.1 of the ALS, with a perfect agreement with previ- 
ous results [lH. This indicates that the crystalline order 
remains coherent along this direction, i.e. perpendicular 
to graphene layers, over a length scale larger than the 
probing depth of ARPES (order of 10 A). 

We now discuss additional features at 2.9 (ei), 4.3 (e2) 
and 7.8 eV (es) characterized by no dispersion at all. 
These features can already be observed in Figure l(b,c) 
taken at 60 eV photon energy, appearing as sharp ex- 
tended horizontal lines at the same energies. Figure 3 
shows a detailed view of these same features observed at 



another photon energy of 55 eV. In panel a, we show the 
first derivative of the ARPES intensity map with respect 
to energy. The first derivative enhances rising or falling 
edges in the data, and thus is particularly useful for de- 
tecting non-dispersive peaks and edges. The energies for 
these non-dispersive features are marked by arrows on 
the right of this panel. In panel b, we show EDCs at T 
point (thin line) and kn= 0.4 A^ 1 (thick line). In these 
panels, in addition to highly dispersive features that we 
already discussed, one can clearly see the non-dispersive 
features, i.e. peaks at energies ei and e3 and edge at en- 
ergy &2- The non-dispersive nature of these features can 
be checked more explicitly using EDCs. For example, in 
panel c, a stack of EDCs over an extended momentum 
range clearly show the non-dispersive nature of peak at 
ei and edge at Such an analysis shows that the fea- 
tures at ei, &i and e3 are non-dispersive within « 200 
meV. 

Similar non-dispersive features have been noted and 
attributed to surface states or emission from isolated 
atoms II 'J. However, it is important to point out that in 
our data there is a strong connection between these non- 
dispersive features and the highly dispersive features dis- 
cussed above. Indeed, one can easily see from Figure 1(c) 
and Figure 3(a) that the energies of these non-dispersive 
features correspond to the band energy extrema. Namely, 
Ci corresponds to the saddle point of it band, e2 the de- 
generate maxima of a\ and <Ji bands, and e3 the saddle 
point of the G\ band. This strongly suggests that they 
are related to the bulk band structure rather than surface 
state. We suggest that these non-dispersive features can 
be explained by loss of electron momentum information 
in two possible ways, both of which are associated with 
the large density of states at the band energy extrema 
|13|. One possibility is indirect transition, i.e. non-k- 
conserving transition in photoemission process |lll Il4| , 
for example, phonon assisted transition. Another possi- 
bility is elastic scattering of electrons in either the ini- 
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FIG. 3: (a) First derivative of raw data taken at 55 eV photon 
energy, (b) EDCs at T point (thin line) and k||= 0.4 A' 1 
(thick line). The positions of these two cuts are shown as 
dotted lines in a. (c) EDCs for the momentum range indicated 
by a horizontal black line (from ki to ki2), marked in panel 
a. 

tial state or the final state by inhomogeneity or disorder. 
In both these scenarios, the observed non-dispersive fea- 
tures correspond to density of states of graphite. This 
explanation is consistent with the EDC line shape ob- 
served in panels (b,c). That is, the non-dispersive fea- 
ture at e2 is sharp at low binding energy side because 
&2 is the energy maximum of <j\ and a 2 bands, whereas 
the non-dispersive features at ei and e3 are broad on 



both sides, because ei and e3 are saddle point energies. 
We note that similar non-dispersive features are also ob- 
served in single crystalline graphite 0, 0] , which sug- 
gests that point defects definitely play an important role, 
in addition to the possible role by the azimuthal disor- 
der we noted above. Interestingly, similar non-dispersive 
features are also observed in NaMoeOi7 ^(| and high T c 
superconductor 17]. Future studies of these features may 
be interesting, in view of the general lack of understand- 
ing of the role of inhomogeneity or disorder in ARPES, 
and the importance of this topic in complex materials. 

In conclusion, the data presented here provide a novel 
example of the coexistence of sharp quasiparticle disper- 
sions with disorder features. Also, several non-dispersive 
features are identified at band energy extrema. Our 
study brings up an interesting possibility that layered 
crystals with azimuthal disorder can be studied with 
ARPES to obtain useful information. 
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